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Abstract 

A theoretical investigation is presented on the properties of charge transfer excitons at vertically coupled 
semiconductor quantum dots in the presence of electric and magnetic fields directed along the growth axis. Such 
excitons should have two interesting characteristics: an extremely long lifetime and a permanent dipole moment. 
We show that wave functions and the low-lying energies of charge transfer exciton can be found exactly for a 
special morphology of quantum dots that provides a parabolic confinement inside the layers. To take into account 
a difference between confinement potentials of an actual structure and of our exactly solvable model, we use the 
Galerkin method. The density of energy states is calculated for different InAs/GaAs quantum dots' dimensions, the 
separation between layers, and the strength of the electric and magnetic fields. A possibility of a formation of a 
giant dipolar momentum under external electric field is predicted. 

Keywords: Magnetoexciton, Vertically coupled quantum dots, Giant dipolar momentum, Galerkin method, Density 
of energy states. 



Background 

Over the last decades, small semiconductor systems with 
discrete energy spectra, known as quantum dots (QDs) 
and that are analogs of atoms, have fired the imagination 
of researchers in many fields of physics [1-3]. Unlike in 
atomic systems, a variety of geometries and configurations 
of the charge state in these artificial atoms are possible. 
These particular features make possible to consider QDs 
as building blocks for the fabrication of more complex 
structures, such as solid-state artificial molecules, case in 
which coupled QDs act similar to coupled atoms in a nat- 
ural molecule [4-6]. Although a diverse range of technolo- 
gies have been implemented to fabricate QDs, in the case 
of artificial molecules, there has been a growing interest in 
spontaneous formation techniques by utilizing self- 
assembling phenomena on crystal surfaces. One of the 
most interesting manifestations of this phenomenon is the 
process of vertical self-alignment of the stacked self- 
assembled quantum dots (SAQDs) [7-9]. These wonderful 
structures composed of two or more vertically stacked 
SAQDs have the advantage of possessing different 
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morphologies such as disks, pyramids, rings or lenses with 
very few imperfections. Also, they are, in general, thin 
layers and have, for the most part, a small height-to-base 
aspect ratio, which is an significant advantage that allows 
us, on one hand, to modify essentially the energy spectrum 
of the particles confined within the heterostructure, mak- 
ing them more stable and, on the other hand, to use sim- 
ple theoretical models. 

Currently, there is significant interest in understanding 
the role of the quantum tunneling of charge carriers be- 
tween vertically coupled QDs, driven not only by a fun- 
damental nature of this phenomenon but also by their 
potential applications. Particularly, many efforts have 
been focused on the theoretical study of the simplest 
configuration of a QD molecule, namely, a pair of QDs 
coupled by tunneling [10-13]. In part, interest in such 
structure arises from its application as a possible gate in 
a quantum processor required to entangle different 
states of an electron-hole pair created optically [14-16]. 
Different exciton states can be disentangled by prevent- 
ing the tunneling through the application of an electric 
field along the growth direction. Formed in this way, 
one of the untangled states, charge-transfer exciton has 
two important characteristics: an extremely long lifetime 
and a permanent dipole moment [17,18]. Additionally, 
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its optical properties can be controlled by means of an 
external magnetic field. 

In this work, we consider heterostructures consisting 
of two vertically aligned hill-shaped InAs/GaAs SAQDs. 
The QDs' morphology has been modeled using a special 
shape which provides an almost parabolic confinement, 
allowing us to perform a relatively simple calculation of 
the exciton spectrum. In our model, a two-dot molecule 
with a single captured electron-hole pair can remain in 
one of two possible configurations with different dipole 
moments. In the first case, when the electron and the 
hole are located at the same dot (on-site exciton), the di- 
pole moment is small, while in the second case, as the 
particles are situated at different dots (charge transfer 
exciton), the dipole moment can be very large. In order 
to illustrate how the electric and magnetic fields applied 
along to the heterostructure growth direction can facili- 
tate or block a transition between two possible carriers 
configurations and in this way control electro-optical 
properties of such structures, we have calculated their 
densities of states and the averaged values of the dipole 
moment for different temperatures. 

Methods 

Theoretical model 

In the case of the charge transfer exciton, we consider a 
model of two vertically coupled InAs QDs in the form of 
the axially symmetrical thin layers with the electron 
located at the lower dot of the radius R e and with the 
thickness at the top W e and the hole located at the 
upper dot with the corresponding parameters R h and 
Wh- In the case of the on-site exciton, both carriers are 
located at the lower QD. Both layers are considered to 
be imbedded inside a matrix of the material GaAs. A 
schematic representation of this system is showed in the 
Figure 1. Here and in what follows, the variables labeled 
by indices e and h are referred to the electron and to the 
hole, respectively. The separation between dots along 
the z-axis is denoted by d. In order to obtain results that 
allows us to analyze qualitatively a transformation of the 
properties of the structure under external electric and 
magnetic fields, we adopt a model of the axially symmet- 
rical hill-shaped QD with infinite-barrier confinement 
and a special shape in which the profile of the dot, given 
by the dependencies of the thickness of the layers w p) 
p = e, h on the distance from the axis p p , p = e, h, are 
defined as follows: 



1 + {Pe/Re 



= ; Wh(/? h ) 



i + (Ph/^) 2 



d+W 



= d+W h /jl+(p/R h f 




p 



Figure 1 Schematic illustration of the cross-section profile of 
vertically coupled quantum dots in the radial direction. Here 
it and ?h are position vectors of the electron and the hole, 
respectively. 



As we show below, such profile provides in-plane 
parabolic confinement for which the mathematical treat- 
ment is significantly easier. For the same reason, we 
adopt a model with infinite barrier confinement, assum- 
ing that dielectric mismatch is much smaller than the 
mismatches of the conduction and the valence bands; 
therefore, the probability for the self-tunneling of the 
particles between QDs is depreciable in comparison with 
the tunneling provided by the external electric field. In 
this way, we assume that confinement potential for the 
electron and for the hole is equal to zero inside the QDs 
and to infinity otherwise, defined in cylindrical coordi- 
nates as: 



V e (r e ) 



Vh(r h ) 



0 if 0 < z e < w e (p e ) 
oc otherwise ' 

0 if d < Zh < d + Wh (p h ) 
oc otherwise 



(2) 



(1) 



Here, three-dimensional position vectors of the elec- 
tron and the hole in the cylindrical coordinates are v p = 
{p p , d p) z p ); p = e,h. 

The values of the physical parameters pertaining to InAs 
used in our calculations are dielectric constant s = 15.2, 
the effective masses in the InAs material layer for the elec- 
tron m e = 0.04ra 0 and for hole m h = 0.34m 0 , the 
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conduction and the valence bands offsets in junctions are 
Vbe = 450 meV and V 0 h = 316 meV, respectively [19]. 

As the quantum dots and the exciton sizes under con- 
sideration are much larger than the unit cell of the ma- 
terial, the effective-mass approximation is a suitable 
approach; therefore the resulting model Hamiltonian of 
the electron-hole pair in the presence of uniform mag- 
netic and electric fields oriented along the z-axis, per- 
pendicular to the plane of QDs, can be written as: 



2^T ( Pp ~^ Ap ) 2 + V p( r p)^P Fz P 



e 

zr e h 



(3) 



p—e^h 



Here, r eh = yj (z e - z h f + r 2 e +r 2 - 2r e r h cos(d e - 8 h ) is 

the electron-hole separation; m p and p p are the effective 
masses and the momentum vectors of the particles, re- 
spectively; and the parameter q p = ± e gives their charges. 
Choosing the gauge for electron and hole vector potentials 
as A e = (B x r e ) / 2; A h = (B x r h )/2, where B is the mag- 
netic field (assumed to be uniform here), the Hamiltonian 
(3) can be reduced to the following dimensionless form: 

H = Hqc + /foh — 2/ r e h; 

H 0p = -rj p [ - A p + y 2 p 2 p /^ - iyspd/ddp+azp^ + V p (r p ); 
p = e,h] rj h = ^/m h ; // e =^/w e ; s h = +1; s e = -1 

(4) 

In these equations H 0e and H oh represent the Hamilto- 
nians of the unbound electron and hole respectively, con- 
fined inside their heterostructures. The following units are 
used in the dimensionless Hamiltonian (4), the exciton ef- 
fective Bohr radius a§ = h 2 e/ft e 2 as the unit of length, the 
effective Rydberg Ry * = e 2 /2s Uq - ft 2 /2^ a§ 2 as the energy 
unit, and y = ehB/2ft c Ry * and a = edoflRy * as the units of 
the magnetic and electric field strengths respectively, with 
{i = m e m h /(m e + m h ) being the reduced mass. The para- 
meters rj e) rj h satisfy the relation rj e + rj h = 1; in our calcula- 
tions, we assume that rj e > rj h . 

Taking into account that typically, in actual self- 
assembled QDs, the height is much smaller than the lat- 
eral dimensions; one can take an advantage of the adia- 
batic approximation, which allows us to exclude 
temporarily the rapid particle motions along z-axis from 
consideration and, in this way, to reduce the dimensional- 
ity of the initial three-dimensional problem. Following the 
adiabatic procedure, described in reference [12], one can 
obtain the two-dimensional effective Hamiltonian, which 
describes only the in-plane particle motion of the form: 



H (2D) =H m + 



j{2D) 
Oh 



rr{2D) _ 
M 0p ~ 

p = e,h 



a {2D) 



2 2 

Y Pp 



+ IP. 

■ iysp 



Phi 



x{z P ) 



(5) 



Here, (z p ), p- e,h are the mean values of z-coordinate 
of the corresponding particle. Finally, for the selected 
profile (1) the Hamiltonian (5) describes two particles in 
2D quantum dot with parabolic confinement: 



p=e,h 

2 



w 2 

p 



■Phi 2 



4tt 2 

W 2 R 2 

p p 



+ r 



(6) 



Finding the eigenfunctions of the Hamiltonian (6) allows 
us to introduce relative and center of mass coordinates r = 
r e - r h and R = rj h r e + tj e r h , respectively. Thus, the 
Hamiltonian (6) is reduced to: 

H = E 0 +H R +H r + U- 

p=e}i 

h r = -n e n h Al D) + (n e <4 + n^l)R 2 /^ 



dd 4 



(7) 



Here, E 0 is the background energy; second and third 
terms H R and H r describe the center of mass and relative 
motions, respectively, and they correspond to the Hamilto- 
nians for two independent central force problems; while the 
last term U presents a perturbation. The Hamiltonian H R 
coincides with one of the circular oscillator, whose eigenva- 
lues can be found exactly and its eigenfunctions can be 
expressed in terms of the generalized Laguerre polynomials. 
It is seen that the Hamiltonian (7) becomes completely sep- 
arable for a particular case when the external magnetic field 
is zero (y = 0), and the geometric para- 

meters of QDs satisfy the following condition: 



(8) 



In this case, 11=0, and eigenvalues of the Hamiltonian 
(7) are a sum of the background energy E 0 , the center 
of mass energy E R defined by two quantum numbers, 
radial N and angular M, and the relative energy E r de- 
pending on four quantum numbers, radial n and angular 
m (N, n = 0, 1, . . .; M, m = 0, ± 1, . . .): 
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E = E 0 + E R (N,M) + E r (n, m); 



£ R (7V,M) = (|M| + 2AT + 1) Jrj^^l + %^h) 



(9) 



The corresponding wave functions are: 

A= \\j {le°>l + f lh CO h)/ f le'1k> 

(10) 

where 4f are generalized Laguerre polynomials and 
< ^ ) w,m( r ) is the radial part of the wave function describing 
the relative coordinates evolution, which is a solution of 
the following ordinary differential equation: 

1 d d0 nm (r) 



r dr dr 

2 

+ 



- +rm + l (n 3 e c>t + nlc>l) r -^f=i 

= E r (n,m)®„, m (r). 



t>n,m(f) 

(11) 



Once Equation (11) is solved and the set of wave func- 
tions (10) is found then it can be used as the basis to cal- 
culate the energy corrections due to the presence of the 
perturbation U in the Hamiltonian (7) in the framework 
of the so-called exact diagonalization or Galerkin method. 

Results and discussion 

In our numerical work, we solve Equation (11) by using 
the trigonometric sweep method [20] initially for d = 0 
and later for d * 0. The eigenfunctions and eigenvalues 
found in the first calculation were used for calculating 
the energy levels E^\ k = 1, 2, 3. . . of the on-site exciton 
(the electron and the hole are mainly situated at the 
same QD), while the results of the second calculation 
were used to find the energy levels E%\ k = 1, 2, 3. . . of 
the charge transfer exciton (the electron and the hole are 
mainly situated at different QDs). 

Once the energies of the two possible exciton config- 
urations E^ and E%\ k = 1, 2, 3 ... are calculated, then 
corresponding curves of the density of energy states /? (0) 
(E) and p {t) (E) can be found by using the following rela- 
tions: 



P (0) (£) = E/( £ - £ f) ; 

k 

f(x) = exp(-x 2 /2s 2 )/sV2n. 



(12) 



Here, the parameter s is a natural width of the individual 
spectral line of the Gaussian shape. In Figure 2, we present 



an example of calculations performed for densities of state 
of the on-site and charge transfer excitons confined in ver- 
tically coupled QDs of radii R e = R h = 80 nm, of height 
W e = Wh = 8 nm and with separation between them 
d = 40 nm. It is seen that left-side energy threshold of 
the curve for on-site exciton is lower than for the 
charge transfer exciton. It is due to the fact that the 
attraction energy between the electron and the hole 
in the case when the particles are located at the same 
QD is higher. 

As for the on-site exciton, the projection of the dipole 
moment over symmetry axis is zero; the external electric 
field does not change nor its energy levels nor its density 
of state p^ 0 \E). On the other hand, the energies of the 
on-site exciton are lower than those of charge transfer 
exciton; therefore, the dipole moment of this structure 
with captured exciton in the ground states is almost 
zero. However, the charge distribution for the ground 
state can be changed drastically under external electric 
field F applied along the symmetry axis which can pro- 
vide a reordering of the energy levels due to a lowering 
of the energy of charge transfer exciton in a value about 
eFd. To verify a validity of this affirmation, we calculate 
the averaged dipole momentum (p) of the vertically 
coupled QDs with a single captured exciton at the 
temperature T by using the following relation: 



ed J p^(E - eFd)exp(-E/kT)dE 
J \p(*) (E - eFd) + p(°) (£)] exp(-E/kT)dE 



(13) 



In Figure 3, we present the calculation results of the 
dipole moment as function of the electric field strength 
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Figure 2 Density of states for exciton confined in vertically 
coupled QDs. The red line corresponds to on site exciton, and the 
blue line corresponds to charge transfer exciton. 
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Figure 3 Electric dipole moment of exciton induced by external 
electric field. Dependence of the dipole moment of the exciton 
captured by vertically coupled QDs on the strength of the external 
electric field applied along the symmetry axis. 



for four different temperatures for the exciton captured 
by the coupled QDs of radii 50 nm, thickness 4 nm, 
and the separation between them 40 nm in the presence 
of the magnetic field y = 3. It is seen for very low tem- 
peratures that the dipole moment increases drastically as 
the external electric field increases and achieves a critical 
value about F c « 7 kV/cm. It is due to the fact that under 
increasing external electric field, all energy levels of the 
charge transfer exciton descend, displacing the threshold 
of the density of states toward the correspondent value of 
on-site exciton. When electric field reaches a critical value, 
the positions of two thresholds are interchanged, while the 
ground state of the exciton suffers a transformation from a 
configuration corresponding to on-site exciton up to a 
configuration of the charge transfer exciton. Such trans- 
formation is accompanied by a drastic growth of the dipole 
moment. In our calculations, we find that the gap AE 



between thresholds of the densities of states of both ex- 
citon types is related to the critical value of the electric field 
F c as AE = eF c d. In this way, values of both F c and AE are 
defined by the competition between the structural confine- 
ment and the electron-hole interaction. 

One can see other peculiarity in Figure 3, the curves 
for different temperatures all have intersections at the 
same point. In other words, when the electric field F 
becomes equal to a critical value F ci the value of the di- 
pole moment is the same for different temperatures. 
One can explain this result, taking into account 
that densities of states both for the on-site exciton and 
for the charge transfer exciton close to their thresholds 
are almost linear. Therefore, when F = F c , the densities 
of states of two types of excitons satisfy the relation 
p^\E - eF c d) = ap^°\E). Here, a is the ratio of the 
slopes of the linear parts of curves close to their 
thresholds. Substituting this relation to Equation (13), 
one can obtain (p) = eda I (1 + a), i.e., dipole moment 
does not depend on temperature. 

An additional possibility to control the properties of the 
structure with captured exciton offers the application of 
the external magnetic field. As the diamagnetic confine- 
ment for the charge transfer exciton is stronger than for 
the on-site exciton, the gap between their ground states 
energies is changed slightly under external magnetic field, 
displacing their energies thresholds in different degrees. 
But more notable is the change which suffers the slope of 
the curve of the density of states of the charge transfer ex- 
citon under external magnetic field. It is clearly seen from 
Figure 4 where we present densities of states calculated 
for charge transfer exciton for two different values of the 
magnetic field. 

Conclusions 

In conclusion, our results demonstrate that electro-optical 
properties of an exciton confined in two vertically coupled 
quantum dots can be changed remarkably by electric and 
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Figure 4 Density of states corresponding to charge transfer excitons. The excitons are confined in vertically coupled QDs for two different 
values of the magnetic field (a, b). The dashed line corresponds to the case without perturbation and the solid line corresponds to the case with 
corrections given by the perturbation. 
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magnetic fields applied along the growth direction. Espe- 
cially, we find that through the electric field, a strong di- 
pole moment can be induced; this is due to the tunneling 
of charge carriers across the potential barrier between 
dots, which leads to a charge redistribution of electron- 
hole pair in the structure, passing from a configuration: 
on-site exciton (the electron and the hole are mainly situ- 
ated at the same QD) to a configuration charge transfer 
exciton (the electron and the hole are mainly situated at 
different QDs). 
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